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Final Report 


ON THE APPLICABILITY OF lOTEGRATED CIRCUIT 
TECHNOLOGY TO GENERAL AVIATION ORIENTATION 
ESTIMATION 


Chapter I 


A. INTRODUCTION 

Current Sensors for General Aviation 

Many of the panel instruments used in general aviation have remained 
unchanged for over half a century. They work and pilots are familiar with 
their use. Considering the retraining required to familiarize nearly a 
million pilots with new displays, it is important that significant valve 
is established before a change is seriously considered. Criteria include 
ease of interpretation, reliability, weight, cost, power, and repairability. 

Advances in Semiconductors/Sensors 

The advances in analog and digital integrated circuits is well known. 
There has been a similar advance in sensor technology which has resulted 
in part from the technology developed to manufacture semiconductor devices. 
Pressure transducers and more recently, accelerometers, have been devel- 
oped which are extremely small and low in cost. The continuation of this 
development is expected. Devices of this type make possible an integrated 
package of sensors which would change the single function sensor and display 
module format of the past and replace them with a single sensing package 
which couid be used to drive a variety of displays and/or autopilot func- 
tions, 

_3. The Estimator 

The other development which can change the approach to orientation 
sensing is in the field of information processing. Many times there is 
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coupling between states needed for display or control which can be modelled 
mathematically. Instead of requiring a measurement for every state, a 
few measurements used in conjunction with the appropriate models can be 
used to generate all the desired states. The formality of this process 
IS extensively documented under the names^ of Kalman Filtering, Estimation 
Theory, or Observer Theory — according to the nature of the problem. 

These theories make it possible to incorporate linear and angular accel- 
erometers, magnetometers, and air data in an efficient way to estimate 
angular and linear velocities, heading and orientation states. Hence the 
information processing makes it possible to incorporate the new sensor 
technology and .to avoid duplication in sensors to > the greatest extent 
possible. It also provides a framework for comparing different models, 
e.g,, kinematic equations vs dynamic equations, or a combination of some 
of each. 

Scope of the Research 

This research has examined some early discussion [Ref, 1] and added 
kinematic equations for comparison. An instrument survey has been per- 
formed to establish the present state of the art in linear and angular 
accelerometers, pressure transducers, and magnetometers. Gyros have not 
been included as they are reasonably well known and one goal of the re- 
search has been to establish if they could be omitted and still obtain 
an acceptable system, A very preliminary evaluation has been done of the 
computers available for data evaluation and estimator mechanization. 

The report develops the theory used, documents the mathematical model 
of a light twin aircraft employed in the evaluation, presents the results 
of the sensor survey and discusses the results' of the design studies. 

In addition to personnel at Stanford: R.A. Van Patten, R. Clappier, 

Russ Hacker, and D. B. DeBra, a subcontract with Stanford Research Insti- 
tute enabled us to benefit from the expert assistance of Dr. M. G. Tashker. 
Our work was coordinated through Dr. Denary at NASA Ames with similar work 
being done by Dr, J. A, Sorensen as an NRG Fellow. This report contains 
portions of Dr. Sorensen’s work. 
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The work was closely integrated; the longitudinal equations and 
estimator represent Dr, Sorensen's work, though there are contributions 
throughout from all parties concerned. 


- 3 - 



II THEORY 


A. GENERAL FORMULATION 

For a general system written in state-variable notation 

X = Fx + Gu + Fw (2.1) 

with control vector u and disturbance vector w and measurements 

z = Hj^x + V (2.2) 

with measurement noise v, an estimator is formed as 

» 

X = Fx + Gu + K (z-z) (2.3) 

For systems where measurements are made of the derivatives of the states, 
the measurement can be written 

z = + v (2.4) 

If the derivative of the state is replaced by substitution from the state 
equation then 

z = (Hj^ + H 2 F)x + H^Gu + H^Fwl v (2.5) 

and the measurement equation contains noise that is correlated with the 
noise in the state equation. Since optimal control techniques require 
that the measurement and plant noise be uncorrelated, a new technique is 
required. 

If H and V are defined as 

H = + H^F (2.6) 
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V - H2rw + V 


(2.7) 


then the measurement may be written 

z = Hx + H^Gu + V 

i 

and the state equation may have added to it a quantity that is 
identically equal to 0, multiplied by a constant L 

X = Fx + Gu + rw + LCz-Hx-H^Gu-v) (2.9) 

which, with the following definition 


F = F-LH 
G = (I-LH2)G 
w = rw-Lv 


may be rewritten 


X = Fx + Gu + Lz + w 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


(2.13) 


If the covariance of the state noise w is Q, and that of the measurement 
noise v is R, then the covariances of the equivalent noises w and v 
are 


Q = E(ww^) = (I-LH 2 )rQr^(I-LH 2 )^ + LRL^ (2.14) 
R = E(w'^) = R+H2rqr^H2 (2.15) 


and L may be chosen to decorrelate the equivalent state and measurement 
— ~T 

noise, E(wv ) = 0, 


— 1 T T 

L = R rqr 


(2.16) 
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An estimator can be built for this system since Lz is known 


A A A 

X = Fx + Gu + Lz + K(z-z) (2.17) 

where 

A A 

Z = Hx + H^Gu (2.18) 

This design can be performed by root square locus, eigenvalue decomposi- 
tion, or pole placement techniques. After the value of K has been 
calculated using the equation above, the estimator can be mechanized 
by resubstituting the definition of F and G which yields 

A A A 

X = Fx + Gu + (L+K)(2-2) (2.19) 

where K has been found from the reformulated problem, and L is given 
in equation (2.16). Equation (2.19) is preferable to equation (2.17) 
because it is not driven by the measurement and because it will have the 
same behavior as the original plant, which is already familiar to the 
design engineer. The gain. (L + K) is the particular value of K which 
minimizes the expected error in the estimator with measurements of state' 
derivatives. 

If the system is designed by pole placement or other methods not 
associated with optimal control, the estimator is mechanized as in 
equation (2.19) with gain chosen to yield the desired performance. The 
following derivation of the estimator error equations holds regardless 
of the method used to determine the estimator gain. In this case, the 

A 

estimator is mechanized with estimated values of the system matrix F, and 

A 

the control distribution matrix G. In addition, it is assumed that there 

A 

are inaccuracies in actual measurements, giving rise to H. Thus the 
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estimator is 


X = Fx + Gu + K(z-z) (2.20) 

A A/ N A A 

z = Hx + H 2 GU (2.21) 

Substituting values of z and z into (2.20) yields 

A A A A A 

X = KHx + [F - KH]x + EG" + KH 2 G]u + KH^fw + Kv (2.22) 
with the definitions 


G" = (I-KH 2 )G (2.23a) 

G" = (I-KH2)G (2.23b) 


Using equations (2.1), (2.22), and the definition of the error in the 
state estimate x=x-x, the error equation may be written 


X = [’(F-KH) - (F-KH)]x + (F-KH)x 

A 

+ (G" - G'')u + Kv - r"w 


(2.24) 


where 

r" = (I~KH2)r (2.25) 

Equation (2.24) shows the dependence of the estimator error on in- 
accuracies in the knowledge- of the plant and measurement. Note that 
in each case, all. (”~) and (") quantities (e.g., H, G") reduce to their 
standard definitions when H 2 = 0 ; l.e., the measurement is not a function 
of the derivatiaves of the states. The effects of the errors in G and 
H will be discussed in later sections as specific mechanizations are 
presented. 



• B. TJKDISTTJRBABILITY 


Breza and Bryson [2] have defined an undisturbable mode as one that 
is not controllable by the process noise, although it may be observable 
through the measurement. Using a Kalman filter such a mode would have a zero 
steady state filter gain indicating that the steady state filter pays no atten- 
tion to the incoming measurement data and hence does not correct the initial esti- 
mate error of the mode. Undisturbable modes can be caused by pure inte- 
grations resulting either from modeling constant disturbances or from 
kinematic relationships in the state model. For example, the lateral 
equations of an aircraft with measurements of the roll and heading angles 
and process noise being a random lateral wind has an undisturbable head- 
ing mode. The filter thus has a neutrally stable eigenvalue which would 
cause the estimate of the heading mode to diverge. 

One solution to this problem is to partition the system into un- 
disturbable and disturbable modes, and then to design an optimal filter 
for the disturbable modes, and an observer with arbitrary dynamics for 
the undisturbable modes. A second method involves introducing aritifi- 
cial process noise into the undisturbable modes. This second method was 
used in this work. The techniques of the previous section were used to 
determine the L matrix producing the reformulated .problem. The w bar 
matrix was then modified so that the random lateral wind drove the state 
equations for roll and heading. 

C. POLE ASSIGI^MENT 

In certain cases, optimal techniques produce estimators with rela- 
tively unpleasant characteristics. Most often, these are manifested by 
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unacceptably long time constants. The engineer has generally produced 
such a design by plugging estimated state and measurement noises into a 
design program, and is naturally unhappy with its consequences. He is 
faced with two choices : either to change the ratio of noises to produce 

a more desirable response, or to abandon optimal control entirely in 
favor of arbitrary pole placement. The former solution is made more diffi- 
cult by the inapplicability of root square locus techniques for systems 
with more than one measurement. The relationship between noise ratios and 
eigenvalue movement is often not obvious. The pole placement technique, 
however, generally presents the designer with an overdetermined set of 
equations for the eigenvalues when there is more than one measurement. 
Estimator gains may be chosen by "closing the loop" with root loci or with 
other computer programs. Pole placement was used in the latter part of 
this work after being faced with slow response of certain variables. 
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Ill AIRCRAFT MODEL 


The aircraft used for this study was the Piper PA-30. The equations 
of motion are those for straight and level flight. 

Lateral Aircraft Model 

The lateral equations of motion in body axes are [Systems Tech. Inc, 176-1, 
page C-3]: 



with the coefficients shown in Table 1. Note that all the aerodynamic 
terms are functions of density and velocity, with those in the control 
distribution matrix being a function of the square of velocity. The 
terms and are the components of the velocity vector along the body 
X and z axes respectively, stability axes not being assumed. 



Table 1 


LATERAL AERODYNAMIC COEFFICIENTS OF THE PA-30 



.3PS.0DUCIBIL1TY OF THE 
-HM-lNAIi PAGE IS POOR 







(3.2) 


The following definitions for the dimensional derivatives appearing in (3.2) 

are used. 

« 
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(3.3a) 
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(3.3j) 
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(3.3k) 

(3.32) 


(■3.'3m) 


(3.3n) 


The dimensionless coefficients which appear in Eqs. (3.3) are given in 
terms of the aircraft stability axes. The coefficients appearing in (3,3f) 
to <3.3n) are thus defined by the following nomenclature: 


C 

zq 


“^Lq 

C . 

za 

- 

'^La 

c 

za 
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c 
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C 
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xu 
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“^DSe 


(3.4) 
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For straight and level flight, the Lift coefficient is computed as 


mg 
QS ■ 


From the value, the trim angle of attack^ can be found from 


(3.5) 


'^L “ ^LO ^ ^La ‘^trim * 


(3.6) 


In this study, for convenience it was assumed that the aircraft body 
axes were aligned with the aircraft stability axes for each steady flight 
condition examined. Stability derivatives were computed for several flight 

t 

conditions; these were subsequently used in the study. Pertinent values 
of the coefficients and derivatives are listed in Table 2. The approach 
flight speed was computed as 


V , = 1.3 V + 0.5 (surface winds) + (reported gusts) 

approach stall 


= 147 + 0.5 (33) + 12.5 (3.7) 


= 176 ft/sec. 

Cruise speed was assumed to be 185 kts (304 ft/sec) . An intermediate speed 
of 165 kts (279 ft/sec) was also used. Altitudes of 0, 1500 ft, and 
5000 ft were assumed at various points in the investigation. 

Other pertinent data for the PA-30 aircraft include; 

= 178 ft^ 

= 5 ft 

2 

= 1900 lb • ft . sec 

= 7.28 

= 0.88 

= 0.33 
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Table 2 


STABILITY DEEIVATIVES OF THE PIPER PA-30 FOR VARIOUS 
FLIGHT CONDITIONS 



P; Ib.s .ft 


30 5000 


0.002377 0.002048 



0.002377 


5000 


0,002048 



0.002377 


Q: lb 


36.81 


79.71 109,83 


94.64 


M : s 

q 


-3,429 -4.684 


-5 . 923 


- 5,104 


-6.516 


-1 -1 I 

M : ft s -0.04660 -0.07668 -0.09697 -0.08355 


- 0.1212 


M ; 
u 


Mr : s 

be 


f-32.62 


-75,95 -100.16 


-88.90 


-120,19 


Z ; ft.s 7.14- 


12.37 


10.71 


15.46 


Z ; s 
I w 


■Cl, 388 


-1.769 


-2.236 -1.927 




-0.3613 -0.2306 -0.2117 


-0.2117 


- 0.201 


Zj. : s 
Oe 


-67.15 


-159.83 -270.28 


-215.64 


-337.85 


X : ft.s 

q 


X : s 

w 


0.10660 0.07099 D. 06516 0.06516 


0.06190' 


X : s 

u 


-0.03061 -0,02191 -0.02667 -0.02298 


X„ • .s 
oe 


Derivative set based on typical flight test uncertainty [4, 5], 
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IV. SURVEY OF INSTRUMENTS 


The survey of linear and angular accelerometer, angular rate sensors, 
magnetometers, absolute and differential pressure transducers, linear 
and angular position transducers , as well as a few temperature sensors 
is essentially completed. Tables presenting relevant data for one or more 
models of each manufacturer are attached as Appendix A. Additional data 
for all models listed^ as well as data on other models made by each manu- 
facturer but-not listed is on file. In some sensor types, silicon tech- 
nology is here, whereas in others, development is still short of the possi- 
ble goal of adequate performance at very reduced cost , Low cost , reliable 
sensors of any type would tisfy the goals of the study. The state of the 
art in reaching that point is variable. Each category of instrument 
warrants some discussion. 

Linear Accelerometers . The available units > surveyed range in 
price from $120 to $595 in small quantity with $50 in 1000 
quantity being the lowest available price. This undamped Entran unit 
is a cantilever beam type with -semiconductor strain gage (piezore- 
sistive) half bridge readout. The very high thermal bias sensitivity 
(5g/l00®F) is one of its disadvantages. This number is more than 
3000 times worse than that of the Teledyne FP 1. Of about 30 t3q3es 
included, only four give data on expected life. At least four com- 
panies making triaxial units were included in the survey; Entran, 

Setra, Humphrey, (the Donner 4384 is single axis), and Donner. The 
Setra 113 triaxial is available at $142 per axis in small quantities. 
Silicon transducers are under development but are marginal at the present 
time. Flight data should be obtained with conventional instruments. 

Angular Accelerometers . Only three companies surveyed make 

angular accelerometers. Of these, one makes units ranged too high 
3 4 2 

(10 to 10 rad/sec ); another (Donner) makes units selling for 
$2500 in quantities of 100 to 500, but the third (Schaevitz) makes 
suitable units selling at about $500 in 100-200 quantities. All 
but the first are listed. 
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Linear accelerometers could be used to measure angular acceler- 
ometers, in principle; but the performance requirements are too exacting 
to make them competitive for the proposed application. For further dis- 
cussion, see App. B. 

This IS the least probable area, in which silicon technology is 
likely to provide a quick solution to high cost at an adequate performance. 

’Angular Rate Sensors . Only three companies have been found making 
rate sensors. The Donner 8160 is presently used in many commercial jet 
aircraft for rate gyro replacement . It is the angular accelerometer with 
electronic integration. The Humphrey devices (manufactured under license 
from Hercules, Inc.) is a true rate sensor using a gas jet passing be- 
tween parallel hot wires. Our data on the B.A.C. Ltd unit is very sketchy 
at present , 

While these rate sensors may be acceptable in performance, the initial 
cost appears too high for general aviation even though it is justified for 
commercial aircraft on the basis of reduced maintenance. At the present 
time, rate gyros may be a better interim solution. 

1 

Magnetometers . Of 15 companies claiming to manufacture magnetomet- 
ers surveyed, only five offer units which are suitable for this appli- 
cation, Mo&t oC the others are for geological sample testing. 

All of the listed units are of the fluxgate type with three of the 
five being triaxial. These three- have a common origin in a rocket flight 
magnetometer application for which NASA Goddard issued contracts with 
a common specification. This accounts for the similarity of these units. 

The lowest priced unit listed (Infinetics) is not satisfactorily cov- 
ered with specifications as its normal use is weapon detection for air- 
port security; however, the manufacturer has indicated that they should 
be capable of our requirements if calibrated. The per axis price is a 
factor of 5 to 6 lower than the triaxial units. 
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Flight test data are needed here more than for any other instrument. 

They have not been used in aircraft as a vector reference before. The 
data available is generally just heading and then that is smoothed with a 
gyro. The more expensive nulled (feedback) instruments have adequate per- 
formance, but may be degraded by aircraft installation problems. It re- 
mains to be seen if the lower cost units will hold calibration adequately. 

New algorithms are to be evaluated, and real data must be used to obtain 
reliable results. 

Absolute Pressure Transducers (altitude). The attached Tables, App. A, 
compares transducers made by 22 companies. Accuracy is given as a percent 
ful scale in various categories. The sea level altitude error was calculated 
.as a root sura squared combination of static error (which includes linearity, 
hysteresis, and repeatability), thermal bias shift for a 4 - 50°F devia- 
tion. In cases where a combined thermal error was given, this was used in 
root sum squared combination with static error as defined above. In 
some cases, the available data was not sufficient, to permit the calcula- 
tion. Items not included in the sea level altitude error calculation were 

long term stability due to insufficient data and initial zero balance 
which would be trimmed out. It was of course assumed that the local 
barometric correction was made by the pilot. Four of the transducers 
listed provide the functionalization required, to. obtain an output iin- 
earilly proportional to altitude. These are the Roseraount 1241A and 
542kl, the Bourns 200-438-1002, and of course the Honeywell HG 280 air 
data computer system. In all other cases, the functionalization is a. 
user responsibility. The Rosemount 542K1 and 542K2 provide both direct 
airspeed (IAS) and dynamic pressure (q^.), in. addition to direct altitude (h) 
the only difference being the airspeed range. Some improvement is possible 
with the low priced sensors by calibrating each sensor against nonlin- 
earity and temperature by least squares fitting. However, the tradeoff 
of this against buying a more accurate sensor needs further study. 

An alternative to the use of an altitude pressure transducer would 
be to increase the resolution of the optical encoder presently used 
with altitude' encoding altimeters by a factor of 10 such that 10 ft 
resolution is available. This approach could lead to a less expensive 
means of obtaining the data. 
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The silicon technology is not here yet for altitude measurement. The 
required improvement in performance is very large, A conventional instru- 
ment must be used unless and until there is a change in the state of the art. 


Differential Pressure Transducers (airspeed) . The Table (App. A) 
compares 21* transducers made by 17 companies. Airspeed ranges for the 
units considered vary; however, in general, most companies will supply 
any range desired. Accuracies are given as static error (defined above), 
thermal bias error, and thermal sensitivity error. Although airspeed 
error is not estimated in the Table it can be calculated approximately 
by combining the sources of pressure error as desired to obtain a com- 
posite pressure error and then using the following relationship: 


error (knots) 


2 

(range knots) PSID% error 10~^ 
airspeed knots ^2 ^ 


For example, a sensor with a 2 PSID (285 knots) range, and 0.5% composite 
pressure error gives an airspeed accuracy at 100 knots of ±2 knots. Three 
of the transducers listed provide the functionalization required to 
obtain an output linearly proportional to airspeed. These are the 
Rosemount 542K2, the Bourns 200-538-1002, and the Honeywell H6 280 air 
data computer. The methods of functionalization, as well as cost and per- 
formance, are dramatically different for these three device-types as 
is illustrated below: 


Device 

Functionalization 

Method 

Accuracy* 
-S0-150°F (knots) 

Cost 

Bourns 

538 (IAS) 
438 (h) 

nonlinear 

pot 

(IAS) ±8k at 100 k 
(h) ±825 ft. S.L. 

$755 

Rosemount 

S42K2 

Analog 

Electronics 

(IAS) ±3k at 100 k 
(h) +47 ft at S.L. 

(q^) 0.5 to 1.5% 

$2,360 

Honeywell 
HG-280 
Digital Air 
Data 
Computer 

Many func- 
tions pro- 
vided. 

Used in DC 10 

(CAS) ± 2k at 100 k 
(h) + 30 ft at S.L. 

(Mach) + 0.005 at 0.35 

$15,000 

to 

$20,000 

• 

at sea level (S.L. ) 
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Low cost silicon sensors should prove acceptable. The airspeed 
display can be nonlinear. There may not be a significant requirement for 
functional fitting because it is the dynamic pressure which is measured 
which is of interest in most cases , 

Linear Displacement Transducers . This Table compares linear displace- 

t 

ment transducers made by 18 companies. The transducers surveyed are of 

various types and cover strokes ranging f rom jf 0.005 inches to 500 inches 

at prices from $32.00 to $700.00, in small quantities. Actuation force 

ranges from 6.5 Ib/in, for some spring loaded units to essentially zero 

for free slug LVDT (linear voltage differential transducer) units. Kavlico 

has applied many t 3 Tjes of LVDT units to commercial and military aircraft 

0 

applications with results like <15 x 10 revenue flight hours without a 
single operational failure. 

Angular Displacement Transducers . Although angular displacement trans- 
ducers were not directly surveyed, many of the companies sent literature 
in angular displacement devices such as potentiometers, LVDTs (linear 
variable differential transformers), encoders, and beam deflection with 
strain gage (linked to shaft rotation). In addition, there are many 
types of resolvers and synchros available for- this application. Exam- 
ples of leading companies manufacturing these components are shown below. 


Angular Transducer Type 

Company 

potentloseter (pot) 

Bourns 

Spectral 

TRWf 

Dale 

LVDT8 

Scbaevitz 


Kaulico 

Encoders ^ 

Astrosystems 

Baldwin 

GUfton/Littoa 

Slnger/Kearfott 

Beam Deflection 

Vest Coast' Research 

Eesolver/Synehro 

Slnger/Kearfott 
Bend lx 

Cllftoa/Lltton 
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A. METHODS FOR TESTING CANDIDATE TRANSDUCERS 

The following methods are planned for testing the candxdate transducers. 
In some cases, the tests are simpler than conventional testing methods. 

Most of the tests will be performed at Stanford, however, facilities for 
some types of transducers are ’not conveniently available on campus and we 
plan to use nearby facilities. 

Linear Acceleration 

We are currently evaluating some accelerometers. We are using two 

testing methods. The first is the standard gravity test signal on a dividing 

-5 

head. A Leitz dividing heas is used which gives 10 rad accuracy and 
hence provides us with two axes of orientation of the instrument with 
respect to gravity and enables us to do some cross coupling and nonlinear 
testing in addition to scale factor and bias. The second testing method 
IS on a shaker which lets us look at nonlinear terms for higher g. We 
also look at harmonic content generated and wave shape distortions for 
large signals. Depending upon the nature of the instrument output, 
more or less calibration>»work will be done on each of these facilities 
depending upon the nature of the information to be sought. 

Angular Accelerometers and Angular Rate Sensors 

We have decided to do^the evaluation of both these types of instru- 
ments on a torsional shaker. Adequately large signals are available on 
the one in our laboratory at Stanford, It has a 0 to 15 deg amplitude, 

at a frequency that can be chosen between 0 and 30 Hz, An optical method 
will be used to determine the peak amplitude. Hence, angular velocities 
an order of magnitude larger than we expected to see in an aircraft, are 
available and considerably higher accelerations than are of interest can 
be generated. Because the angular accelerometers and rate sensors are 
to be used in combination with magnetometers and accelerometers which will 
be the primary source of long terra orientation information, no conventional 
gyro drift test will be performed. 
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Magnetometers 


Testing magnetometers require a large, complicated facility. Stanford 
has no such facility but there is a very good one at NASA Ames that we 
should be able to use. There they are able to measure magnetometer zero 
offset, linearity, and drift; and they routinely run tests on flux- 
gate magnetometers of the type we are considering. 

Air Data Sensors 

Air data sensors, including absolute and differential pressure trans- 
ducers, are usually tested using precision servo controlled pressure 
sources. There are several adequate test facilities in the Bay Area, 
including one at NASA Ames, If necessary, an acceptable system could 
be purchased for $10,000 to $15,000. Another alternative would be to 
use the local instrument shops that serve the general aviation community, 
but their equipment may not have the required accuracy. 

Linear Displacement 

The work table of a sip jig borer, which is calibrated to very high 
precision, is a very convenient work table on which to do displacement 
calibration work. We have had very good success in calibrating LVDTs 
using this setup in the past. - 

Angle Transducers 

Conventional dividing heads should provide adequate accuracy for 

calibrating angle transducers . If interest develops , a Leitz dividing 
—5 

head good to 10 rad is available for especially fine measurements. 
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t B. MICROPROCESSOR INTERFACE CONSIDERATIONS 


The interface between the aircraft' sensors and the microprocessor 
can be either high level d-c analog or ‘bit parallel 'digital. For 
analog signals, an A/D converter controlled by the microprocessor must 
be used. To digitize more than one analog input, a processor controlled 
analog multiplexer must be connected to the A/D converter. Most A/D 
converters require a fairly high level input of 0 ± 1 V to 0 ± 10 V. 
Sensors that produce signals smaller than this will require some signal 
conditioning prior to conversion. Sensor that produce a-c signals 
will probably require demodulation to produce a d-c signal for the 
converter. Actually most A/D converters are fast enough to do the 
demodulation under software control but this often uses a lot of 
processor time and is undesirable unless the processor is lightly 
loaded. In any case, the hardware for d-c or a-c signal conditioning 
IS straightforward except for sensors that produce very small signals. 
Sensors with digital outputs can be connected directly to the micro- 
processor's l/O bus through an appropriate buffer. This buffer gates 
the data on to the 8 bit to 16 bit wide l/o bus under processor control. 

The actual hardware that must be purchased and/or developed de- 
pends on the number of sensors and the type of |j,processor used. 

Some of the converter-producing companies sell complete data acqui- 
sition systems that plug directly into the chassis of some established 
^processors such as the popular Intel 8080. These printed circuit 
boards are completely compatible with the ^processor and require no 
hardware development by the user. A typical system with 16 analog 
inputs and 12 bit digital resolution costs about $700 in single 
quantities and $300 for 100 units. This class of jiprocessor uses 
8 bits words and has limited arithematical capability and may not be 
powerful enough for our application. 
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The recent more powerful uprocessors such as the one used in the 
H-P 9825 and the Data General micro NOVA approach minicomputer sophis- 
tication. Processors such as these used 16 hit words and have power- 
ful arithmetic logic, but compatible plug-in data acquisition systems 
are not yet available. What is available, however, are digital l/O 
boards which provide a well defined, fairly simple digital interface 
to the user. To these boards can be connected A/D converters, multiplexers, 
etc., with a small number of logic integrated circuits. Some converter 
companies produce data acquisition modules that contain the multiplexer, 
the A/D converter, and some control logic which can be connected to 
the digital l/O boards with even less logic. In either case the user 
must do some logic design and is responsible for the compatibility of 
the digital and analog subsystems. Digital iJq boards or modules 
are supplied by the ^.processor manufacturer and cost $200 to $400 in 
single quantities, A typical 16 analog input, 12 bit data acquisition 
module costs $300 in single quantities and about $150 for 1000 
units. 

The amount of hardware design required for such an interface 
depends strongly on the data rate. The above discussion assumed that 
data was gathered in the simple programmed l/O mode. This mode re- 
quires simple hardware but uses a lot of processor time. Data rates 
of 10 to 5000 total samples/sec (e.g., 1 to 500 samples/sec for 
each of 10 sensors) are' reasonable depending on the processor work ♦ 
load. To define a tighter data rate range would require knowledge 
of a specific processor and the software tasks. The program interupt 
mode of data gathering requires somewhat more interface hardware but 
uses less processor time and thus allows higher data rates . The 
Direct Memory Access (DMA) mode requires still mo^e hardware but 
uses very little processor time and allows data rates up to the 
maximum capability of the A/D converter — up to 100,000 samples/sec. 

Most jiprocessors support all three modes of data input. The appro- 
priate mode depends entirely on the application. The interface for 
any of these modes uses standard integrated circuits and the entire 
analog input subsystem is far less complicated then and should be at 
least as reliable as the jiprocessor itself. 
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V 


STATE ESTIMATOR DESIGN FOR IHE LONGITUDINAL MODE^ 


In the longitudinal mode, the two aircraft state variables which are 
not considered to be directly measureable are the pitch attitude 0 , and 
the pitch rate q. The objective of the state estimator is to determine 
estimates of these variables as well as to produce smoothed values of the 
other variables which are directly measured. The smoothing is desired to 
reduce the effects of instrument noise and wind disturbances on the sensor 
signals. 


A. OBSERVER DESIGN 

The linearized aircraft perturbation equations for longitudinal motion, 
including wind and instrument dynamics, are given in Eqs. (3.2). In 
abbreviated form, these equations are represented by the matrix differential 
equation 

X = Fx + G^u + FTlj^ . ( 5 . 1 ) 

Here, 

X = state vector, 
u = control input vector, 

= wind noise vector, 

F = system dynamics matrix, 

G^ = control distribution matrix, 

r = wind distribution matrix. 


The material of this chapter was taken from Ref. 8. 
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The estimator model as 


X = tx + t K(y - y) . 


(5.2) 


Here^ 

X = estimated state vector^ 

u = measured control input vector, 
m 

F = assumed system dynamics matrix, 

G = assumed control distribution matrix, 

and measurement y is 


y = Hx + Tl^ • (5*3) 

H as the output distribution matrix, and indicates measurement noise. 

To obtain an error signal, an estxmate of the measurement y is formed? 

5 ? = hS . (5.4) 

In almost all flight control appli^tions, it is necessary to have 
the estimator produce 0, q, and h. A full state estimator in Aq, Aw, A^, 

Au, plus altitude perturbation Ah is therefore desirable. In addition, 
it is' possibly desirable to estimate winds, (wx, wz), lagged instrument 
measurements, (^^m^ ^m^^ instrument biases, ^^22:^ ^^e^* 

The full ten-state estimator mechanized here used 10 gains with 3 
measurements: 4 gains associated with the vertical accelerometer, 4 with 

the measurement of airspeed, and 2 with altitude measurement. In deter- 
mining the estimator gains, the coupling between the short-period and 
phugoid equations was ignored; the associated gains were selected inde- 
pendently. The coupling was included in the actual estimator mechanization 
however. 

The equations of the short-period mode, phugoid mode, and altitude 
kinematics were combined along with equations for wind components and certain 
biases to produce the ten-state estimator where 
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(5.5) 


(5.6) 


The ten state variables for thxs estimator are: 

= [Aq, Aw, A0, Au, Au^, W^, Ah, Ah J . (5.7)- 

A five-state estimator without wind or bias terms was also tested 
extensively. 

Runs were made to compute and plot the transient state estimate 
errors x (x = x - x) due to modeling errors and estimate initial con- 
dition errors. Modeling errors were simulated based on stability deriva- 
tives in the observer set for a different altitude and airspeed than those 
set for the actual aircraft flight conditions. In one set of cases, 

V of the estimator and the associated derivatives were based on an 
a 

assumed nominal airspeed of 100 kts. , The nominal/ -V (¥„ ) was assumed 

a ®nom 

to be 185 kts. Thus, tihe stability derivatives used to model the actual 
aircraft motion were based on airspeed of 185 kts. The airpseed perturba- 
tion Au which was assumed as the measurement to this system was Au = 
Plots of transient errors A5, and Ah due to errors in 

the observer estimate , initial conditions with these gross modeling errors 
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revealed that although the ten— state observer produces a smaller pitch rate 
error Aj, the pitch angle £0 and attitude errors are smaller for 

the five— state system. The reason for this is that when there are extra 
wind states with incorrect models, the observer takes measurement differ- 
ences and partially converts these to wind estimates. The conclusion is 
that extra states Improve performance only if the modeling is correct. 

For transients due to an elevator input sequence consisting of a 
3° doublet over 0.8 sec followed by a 0.3° step input lasting through 
10 sec, the five-state estimator produces superior performance in 
estimating the vertical and longitudinal airspeed components (Aw, Au) 
and pitch angle W* Pitch rate Aq and altitude Ah errors are 
about the same for both observers. Again, the five-state observer 
proves to be overall more tolerant of modeling errors. 

From results with estimators containing modeling errors, the follow- 
ing general conclusions can be made: 

1. Observers with fewer states are less susceptible to modeling 
errors. Adding states to the observer to increase its accur- 
acy should only be done when a good overall model is known 
and used. 

The pitch angle error transients for the ten-state observer 
with initial estimate errors and elevator inputs are totally 
unacceptable. The pitch angle error transient for the five- 
state observer was excessive only during a gust. Thus, it is 
recommended that the five-state observer be used as a first 
choice flight for test purposes. 

2. The stability and control derivatives are a function of air- 

• « 

speed. It is recommended that the values used for these deriva- 
tives be updated to match the airspeed as closely as is practical. 

The effects of mismodeling were very serious — the stability 
derivatives varying with velocity, and the control derivatives varying 
with the square of velocity. To reduce these effects, the mechanized 
values of the derivatives were changed with airspeed and altitude (all 
derivatives are functions of density) and the following results obtained. 
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' 1. Dividing the flight speed range into finite regimes is an 

effective way to reduce modeling errors'. It provides a simple 
way to implement gain changes. For the PA- 30 ^ dividing the 
span of flight speeds into three segments is adequate. Altitude 
change effects are not important. 

2. The control derivatives are more sensitive to speed changes than 
are the other parameters. It is recommended that these quan- 
tities be updated for every 10 ft/sec change in flight speed. 

3. The effect of uncertainties in the dimensionless stability 
derivatives can be significant, especially on the pitch angle 
error. Careful modeling is necessary, and the effects of 
known uncertainties on the results should be checked. This 
factor is a strong motivation for careful parameter identifi- 
cation of the derivatives. This is most important for the 

control derivatives. Me and Z„ , 

’ Oe oe 


B. KINEMATIC FILTER DESIGN 


For the longitudinal mode of the aircraft, body-mounted transla- 
tional and rotational accelerometers would measure the perturbation 
quantities 


Aa = Au + + g cos 0 A0 + T| ; 

xm o 'ax 

• 

Aa = Aw - V Aq + g sin 0 A9 + T] ) 

zm a ^ ° o ‘az 

Aq = Aq + Tj. . 
m . q 


(5.8) 


Accelerometer measurement errors are indicated by the terms T| ^ ^ 

ax az 

and . Note that these equations contain primarily the kinematic 

terms V , and 0 ; the aerodynamic forces and moments are not explicitly 
a o 

modeled. These equations can be transformed to 
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Au = Aa - , - g cos 0 A0 , 

Aw = Aa + Aq - g sin Q A0 , 
zm o 


(5.9) 



Here, the error terms are neglected. These equations form the basis 
for constructing an alternate state estimator which does not require 
specific aircraft modeling as with the observer approach. It is noted that 
mechanization of this state estimator approach does not have to be based on 
small perturbation' assumptions , as with the previously discussed observer 
approach. However, the small perturbation assumptions and notation are 
maintained to allow linear analysis and direct comparison between the two 
approaches. 

To these equations, the kinematic expressions 


A0 = Aq 

Ah = -cos 0 Aw + Vh A0 + sin 0 Au 
o o o 


(5.10) 


are added. Here, Vh is V cos Q , These are the equations of 
■ o a o 

the state estimator considered here. In (5.9), the measured acceleration 
(Aa , £Si , and ) are treated as inputs (like the elevator position 
inputs for the observer approach). Because (5.9) and (5.10) are made up 
of measured accelerations and kinematic terms, this filter concept is re- 
ferred to as a "kinematic filter." 

In matrix form, the kinematic filter is implemented as follows: 








(5.11) 
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Note that the airspeed measurement Au is not used to update the alti- 

m 

/% 

tude estimate and the altimeter measurement Ah is not used to 

^ ^ ^ m 

update estimates of Aq, A0, or Au. This naturally follows from the 
orthogonality of the equations. Thus, only six gains are required to 
implement this state estimator. The gain K can also he eliminated, 
and this results in the gain selection for the vertical motion being 
decoupled from the forward and pitch motion. 

The obvious advantage of the kinematic filter over the observer 
is that no aerod 3 mamic stability derivatives need to be modeled. Also, 
elevator deflection measurements (6^) aren’t required. The immediate 
disadvantage, from the mechanization point-of-view is that angular 
accelerometer and forward linear accelerometer measurements are required 

Testing this kinematic filter resulted in the following conclusions 

1. The dominant disturbance affecting the estimate was the longi- 

tudinal wind gust. The kinematic filter produces about the 
same error for Au as the observer, slightly larger errors 
for Aq and and substantially larger errors for Aw. 

These error sensitivities are, of course, dependent on the 
filter gains, and a tradeoff exists between steady state error 
and transient response. 

2 , Of the instrument noise errors, only the angular accelerometer 
noise produced a significant effect. Some prefiltering may be 
warranted. 

.3. The only significant modeling error that can occur in the kine- 
matic filter formulation is using an incorrect value for the 
nominal airspeed. • In all cases, errors resulting from an 
incorrect mechanized value of airspeed were only slightly 
greater than if the correct value was used. The kinematic 
filter is much less sensitive to modeling errors than the 
aerodynamic observer. 
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VI LATERAL ESTIMATOR 


The lateral equations of state are shown in -equation (3.1). The 
states are lateral velocity, roll rate, yaw rate, roll angle, and head- 
ing. Currently, steady-state heading is determined using a magnetic 
compass, and short-term heading changes are displayed by a directional 
g 3 rroscope. In IFR flight, the directional g 3 ?ro is reset by the pilot to 
agree with the compass every 10 to 15 minutes in level flight, and before 
each approach to landing. The primary attitude reference is the artifi- 
cial horizon that displays pitch and roll angles. Another gyroscope dis- 
plays rate of turn (yaw rate) . The simplest of all instruments, a damped 
ball in a curved tube, indicates lateral coordination. Ignoring this 
passive Instrument the pilot’s lateral displays are roll angle, heading, 
and yaw rate. 

Aerodynamic Model Observer - 

Quantities that can be sensed for use in the lateral estimator are: 
heading, yaw rate, yaw angular acceleration, roll angle, roll rate, roll 
angular acceleration, lateral velocity, and lateral acceleration. The 
scope of this work requires that angular velocities not be considered, 
they being measurable primarily with gyroscopes. Lateral velocity will also 
be ruled out as it is measurable primarily with lateral angle of attack 
vanes; in addition, the design of an observer would not ordinarily choose 
lateral velocity as a measurement, inasmuch as it does not play a major role 
in the observability of the system. In considering which (if any) of the 
acceleration measurements should be used, the price and complexity of the 
instruments must be taken into account. Linear accelerometers are the least 
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complex and costly, and inasmuch as lateral acceleration is used as a mea- 
sure of coordination, it is a reasonable measurement. The inclusion of roll 
or yaw angular acceleration will depend on which angles measured. If both 
roll and yaw (heading) are measured, both angular velocities are determined, 
and thus, the measurement of the accelerations will just add a high frequency 
component to the observer. If roll is not measured, the measurement of roll 
acceleration will add important information to the observer as roll xrf.ll be 
derived only through the other observed equations. Inas much as yaw must be 
measured to provide the "outer loop" measurement to the entire observer, it 
is reasonable to measure roll rate. Thus, for an observer in which the 
d 3 Uiamlcs of the system are modelled, the observations chosen are heading, 
lateral acceleration, and roll acceleration. 

The equations of such an observer are 
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with three measurements: 


z- = a , the lateral linear acceleration, 

1 y ’ 

Z 2 = ot^j the roll angular acceleration, 
z^ = the magnetic heading. 

Note in equation (6.1) that since measurements 1 and 2 are derivatives of 

A 

states, the control appears in the term K(z“Hx). Such an observer may be 
designed using methods of optimal control by first calculating the decoup- 
ling matrix L of Section II from the estimated measurement and plant noise 
covariances. The equivalent system (F,G,H) may be used in an eigenvalue 
decomposition program such as OPTSYS to find the optimal gain matrix K. 

The actual gain used in the estimator is (K+L). This procedure was -per- 
formed for the PA-30 aircraft for a variety of covariances at an airspeed 
of 176 ft/sec. At this speed, the roots of the aircraft' are: 

' -1 

-3.96 sec , roll subsidence 

-5.49 ± 2.14 sec Dutch roll 

0.00233 sec , spiral divergence 

and a root at the origin representing heading. 

An observer with good transient behavior has roots at: -1.64, -1.02, 

±.7-59, -.140 + .150.sec~\ 

Perhaps as important as behavior of the estimator to transients and 
wind gusts is its behavior when the aircraft receives a control input from 
the ailerons or rudder. This can be considered separately from transient 
behavior due to superposition of responses. The response of the estimator 
error, equation (2.24), to control inputs, is obviously dependent on the 


- 34 - 



knowledge of G''=(I-KH 2 )G as well as the actual control u. It can be seen 
from Table 1 that the elements of G are a function of d 3 mamic pressure; i.e., 
air density and the square of speed. Thus, any error in the knowledge of 
the aircraft's speed with manifest itself in an error in the mechanization 
of G and thus give rise to an effect on the observer by a control input. 
Unless the mechanization updates the aircraft speed constantly and recal- 
culates the mechanized values of G, this error will exist. In addition, 
the terms in the elements of G that are aircraft dependent are typically 

the least understood theoretically of all stability derivatives although 
empirical data is reliable. 

To determine the effect of this mismatch, a 15 second digital simula- 
tion was performed in which the aircraft was subject to a single square- 
wave input in aileron of amplitude 0.1 radian and 10 second period 
(with no forcing term for the last 5 seconds). The aircraft speed is 
176 ft/sec; the estimator is mechanized assuming an aircraft going at 
150 ft/sec. 

Figures 6. la-e show the response of the aircraft along with the esti- 
mator errors. The errors are clearly unacceptable. They can be explained 
by examining equation (2.24). The observer error is a function of the mis- 
match G (which is a function of the square of velocity) multiplied by ob- 
server feedback gain. The term KH 2 (G-G) may be large (N.B. this large 
term exists only in an observer with nonzero H 2 ; i.e., where there is ob- 
servation of state derivatives.) Parenthetically, it can be noted that 
the effect of this term reduces to zero 'if (I-KH 2 )= 0 , a condition that is 
unlikely to be met in a stable system. Work by the experimenter has failed 
to reveal a system with acceptable transient behavior (roots) in which the 
response to control modeling error is small. 
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MTERAL VELOCITY (ft/sec) 


5.5 

' 3.86 




ROLL ANGLE (rad) 





Kinematic Observer 


An alternative mechanization is one in which an accelerometer is used 
as input to a kinematic equation, and these equations are used in place of 
the dynamic state equations. When this is done, the measured value itself 
IS treated as an input, and hence is known just as the control inputs are 
assumed measured and therefore known in the case of the dynamic equations. 
There are certain combinations of linear and angular accelerometers that 
can be used as inputs rather than measurements that will simplify the mech- 
anization and more importantly, reduce sensitivity to modeling errors. 


This set of possible mechanizations assumed that the magnetometer is 
used to generate both heading and roll ‘angles (with knowledge of pitch angle 
from the longitudinal estimator, and the dip angle of the earth's magnetic 
field). In addition, there are a maximum of three measurements; a , lateral 

y 

acceleration, ct , roll angular acceleration, and Ci , yaw angular accelera- 
X y ^ 

tion. Given the assumption that all mechanizations will include the linear 

accelerometer measurement of a due to the lowest instrument cost, there 

y » 

are four possible mechanizations, labeled "Types I through IV." 


It should be noted that extracting both heading and roll information 
from the magnetometer would require either some form of vertical reference 
based on specific force measurements which are filtered to get rid of accel- 
eration (employing a gyro) or a processing scheme that would attempt to re- 
solve ambiguities in heading as the aircraft passed through magnetic north 
or south. The following estimator mechanizations were studied to examine 
the results of using accelerometer signals as inputs and of not using con- 
trol surface deflections in the estimator. 


Type 1 Estimator 


v = -Ur + Wp + (g cos 0)£ + a ' 


p = a 


X 


r = a 


+ K 


$ = p + (tan 0)r 
ij) = (sec 9)r 


^1 


Z2 - Z2 


( 6 . 2 ) 
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where z, = <j> ,, and z„ = i|i ,, and a , a , a are accelerometer 

1 measured 2 ’^measured y x z 

inputs. Note that no control input is required since its effects are mea- 
sured by the accelerometers. In this set of equations, v is not observ- 
able, a drawback for an autopilot, but not for a display. Ignoring the v 
equation, the system can be mechanized as^ two decoupled second-order systems 
in (0,p) with measurement of <ji and input ci , and ('/',r) with measurement ip 
and input The disadvantage here is that two relatively expensive angular 

accelerometers are needed. If use of either one (or both) of the angular 
accelerometers is to be avoided, an equation with aerodynamic terms must be 
substituted for it. The mechanization without the angular accelerometer 
in yaw is labeled Type II. 


Type II Estimator 

Here, aerodynamic terms are used instead of a , thus reducing the in- 

z 

puts to a , a 
y X 




v = -U r + W p'+ (g cos 0)$ '+ a 
o o ® ^ y 


p = a 


r = Nv + Np+Nr+N5, 5a + N~5 
V p r oa 6r r 


$ = p + (tan 9)f 
$ = (sec 0)r 


+ K 


Zi - z^ 


^2 " ^2 




(6.3) 




Note here that the r equation has the same terms as in equation (6.1) 
including control terms. There will be none of the large errors caused by 
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mismodeling of G, there being no state derivative measurements (the terms 
treated as inputs here), in fact, it has proved possible 
to ignore the effect of control entirely, even for control-forced behavior, 
if the gains are chosen appropriately. If the gams chosen are 


K = 


0 48000. 

200. 0 

0 371.^ 

20. 6 

0 29. 


(6.4) 


then the roots of the estimator equation are; -4.08, -12.9 ± 11.2, -10.0 
±10.0 sec While there was no need for the roots to be so fast, this was 
done to simulate the effect of using no control input in the observer. 

Note that these roots were obtained with only five nonzero gains. Figures 
6.2a-e show the results of the same type of simulation done for the aero- 
dynamic estimator.^ There was, however, an assumed difference of 50 ft/ sec 
between the actual speed of the aircraft, and that used to calculate N , N , 

V p 

and N^, the only speed dependent terms in the observer and 

were assumed zero, and thus there would be no need to instrument the 
ailerons and rudder). Note that the only objectionable error is that in 
the lateral velocity, ^^7hlch is approximately 20-30% of the magnitude of the 
velocity itself. Simulations involving a step in rudder yielded equally 
good results for yaw and roll angles and rates, but worse errors in lateral 
velocity. These approached the magnitude of the lateral velocity for the 

Comparable to Figures 6,la-e, 
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LATERAL YELCCITY (ft/sec) 



ROLL ANGLE (rad) 



Fig. 6.2c 


Fig. 6. 2d 






same assumed 50 ft/sec speed error. These errors can be reduced with gain 
scheduling or with the inclusion of the control terms. 


Type III Estimator 

Here, aerodynamic terms are used to replace a 

X * 




v=Ur+Wp+ (g cos 0)^ + a 
o o y 


p=Lv+Lp+Lr+L 6 +L 6 
p r aa rr 


r = a 


2 - - 2 . 


+ K 


(6.5) 


$ = p + (tan 0)r 
\p = (sec 0)r 


2 -, - 2 . 


The Type III estimator, employing the angular accelerometer for yaw 
instead of roll, did not perform as well as the Type II estimator did 
for lateral velocity (without inclusion of control terms) . It did as well 
for the other states. 


Type IV Estimator 

Here aerodynamic terms are used to replace ‘ and a^: 


v = Ur + Wp + (g cos 0) 4) + a 
o o y 
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The Type IV estimator, with only a linear accelerometer, is very close in 
form to the purely aerod 3 naamic estimator, although there is no measure- 
ment of state derivatives. Its simulation response suffered severely 
from lack of control inputs. This is intuitively reasonable since 
there are no angular accelerometers to measure the immediate direct 
effect of the control. Were such a mechanization to be used, both lateral 
control surfaces would have to be instrumented. 

Vector Mechanization 

It IS convenient for conceptualizing and checkout to keep as many 
modes separated as possible. On the other hand, the evaluation of the 
interaction between lateral and longitudinal behavior is best done 
with the integrated equations. For example, pitch influences the magnet- 
ometer readings, and hence heading information needed In the lateral 
equations, A set of combined equations for the orientation is derived 
for any two known spatially fixed noncolinear vectors. Their time 
derivatives in space is zero so in body axes we get 


b 


B = 

^b 

-uj 

X B + 

Kg(B™ - B) 

(6.7) 

b 

/S 

g = 

^b 

‘a 

X g + 

Kg(-f”’-g) + K^bCb"' - B) 

(6.8) 


and Ij 

W = a"* + K „(b"* - B) + K (-f"* - g) (6.9) 

where Kgg and may be zero and b implies differentiation or coordi- 

nates with respect to the body. (Other estimator feedback combinations may 
prove to be more effective,) 

These equations could be simplified slightly by the direct use of 
a measured u from a rate gyro set. Similarly, we can write the vector 
form for the acceleration but using the specific force as an input 

y = X V + + g + - U) + - B) (6.10) 

where again, K may be zero. 

Vj3 
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Finally, the altitude equation using a local frame Jl, for vertical vel- 
ocity component 


^ v\ in ^ 

“ - -'^3^ + V'' - « • 


( 6 . 11 ) 


Although all terms of the direction .cosine matrix are not needed 




known 


from which the necessary terms for 





(6.13) 


can be expressed in terms of and g^. 

Using the same relationship^ the indication of pitch and roll altitude 
will be the direction cosine elements 


pitch 

= -lb . 

3j& 

= - C 

s. 


/\ 


/s 

roll 

= +2b • 

3j6 

= c 


13 

23 


(6.14) 

(6.15) 


which can each be expressed more simply in terms of and g^. 

Since the local frame will be defined as = magnetic north only 

the dip angle need be estimated in flight (or set occasionally) . 


This mechanization has not been simulated in detail. Preliminary 
evaluation suggests that the gams should be appropriate for filter— 

me electronic noise ^ 1 sec whereas the K should be comparable to 

s 

the evection gains for an artificial horizon producing a few hundred 
second time constant and having a cutout for large maneuvers or at least 
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a saturation level for the whole term involving The other gains 

follow a pattern discussed previously in Ch, V. 
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VII 


CONCLUSIONS AND RECOMMENDATIONS 


Modeling aircraft dynamics can give more flexibility in the choice 
of sensors needed to measure adequately the states of an aircraft. The 
models make it possible to infer states with indirect measurement. We 
have shown that an instrument set composed of linear and angular accel- 
erometers, mangetometer, and pressure transducers is adequate. This com- 
plement of instruments should be able to take advantage of the emerging 
state of silicon technology spawned by the semiconductor market which 
promises great reliability and low cost. The use of the aircraft dynamic 
equations in an estimator has been shown to be feasible, but errors are 
introduced which depend on disturbing forces and torques (due to winds) 
and in some axes, it is important to include control inputs. Purely 
kinematic equations are not as sensitive to modeling omissions but there 
are questions of observability. Tims' we have been led to a combination 
of kinematic and dynamic equation relations which minimizes the depend- 
ence on wind and disturbance modeling~and do not require control inputs. 
This is true for both longitudinal and lateral equations. 

In addition, states representing orientation, velocity, and altitude 
in an integrated vector “format have been considered but not evaluated in 
detail. An estimator was designed this way which has the appeal of hav- 
ing a close analogy with an artificial horizon. On the average, the 
specific force is gravity if the velocity is to be bounded. The strong 
dependence on the angular velocity in, these equations suggest that an 
alternate configuration of sensors should include rate gyros. 

Thus it IS concluded that the combination of three accelerometers, 
two angular accelerometers, a three-axis magnetometer, absolute pressure 
and pitot pressure are an adequate set of sensors for the estimation of 
orientation, altitude, and velocity states. Reasonable performance re- 
quirements are required and can be met by existing instruments though not in 
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all cases by the smallest or least expensive as originally hoped. Since 
the new technology is expected to continue to develop, it is recommended 
that the feasibility of this approach be established with available 
sensors on the expectation that reductions in cost and improvements in 
reliability will follow soon and may be available by the time a flight 
evaluation is under way. Silicon technology is adequate for pitot 
pressure now but development is still needed for barometric altitude 
measurement, Silicon technology accelerometers may be acceptable but an 
experimental evaluation must be completed with instrument tests before a 
decision can be reached. Feedback type fluxgate magnetometers are adequate 
but non-feedback types must be tested to determine if their characteris- 
tics are stable enough for online modeling. 

It IS recommended that the candidate instrument package flown -include 
three axes of rate gyros angular accelerometers, linear specific force, 
magnetic field, and the two pressures. In this way extra measurements 
may be included or not in the off-line evaluation. 
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APPENDIX A 


This Appendix includes detailed Tables, mentioned in Chapter 
IV, for the following: 

LINEAR ACCELEROMETERS (4 pages) 

ANGULAR ACCELEROMETERS 

ANGULAR RATE SENSORS (2 pages) 

MAGNETOMETERS 

ABSOLUTE PRESSURE TRANSDUCERS (4 pages) 

DIFFERENTIAL PRESSURE TRANSDUCERS (airspeed) (3 pages) 
LINEAR DISPLACEMENT TRANSDUCERS (3 pages) 



Mfr. &. Life 


Model 


llantje 


Type 






itxvity 

'EntxAn 

EG-D240 

-10 

± lOg 

Piezore- 
sistive j 

C. beam 

10 mV/g 
at lOV 
ex. 

Entran 

EGC- 

240-D5 

±5g 

It 

j 

25 mV/g 
at 15V 
ex. 

Kistlex- 

Morse 

201-020 

±20g 

Piezore- 

1 

sistlve I 

30 mV/g' 
at lOV 
ex. 

Genxsco < 
10^ cycles 

2386 

±5g 

Pot. 

IV/g @ 
lOV ex. 

Genisco 

4 

10 hours 

2388 

±5g 

Variable 

Reluctance 

ac phase 
sens it V, 

C.J. 

Enterprises 

1 

A18 

±5g 

t« 

If 

Timex 

AHI-000 

±2.5g 

Differen- i 
tial i 
xfmr, ! 

1 

tl 

Kulite 

GAIl-813 

-10 

±10g 

semloon-i 

strain 

gage 

12 mV/g 
§ lOV ex 

Teledyne 

FP 1 

±l5g 

force 

rebalance 


G.E. 

1 

LB-5 1 

±20g 

tl 

lOOmV/g 


LINEAR A 


Non- 

linearity 

Thermal 

Sens, 

Shift 

2%FS 

10-16% 

per 

100“F 

1% of 
reading 

+2,5% 

per 

100“F 


Qty Size or 

Price Weight 


5 g per 
100°F 




1-4 @ 
$ 120 ; 
1000 
$50. 


1-9 @ 
$270; 
25-49 
$230. 


1,$160; 
25, $120. 

50 gm 

%$412. 

3”xU" 

150 up 
$337, 

X li" 

1§$563. 

2” X 2 

150 up 
$457. 

X 2^' 



±1% FS 
static 
error 


1%/100°F 0.5g 0,01 

e/s 


1.25X 
— 3 

10 g @ 
20*C 


2% per 0 . 33g 

100°F SO^F to N.S. 5% max. 
180”F 


10 10 ^g/g^ 




vita. Sold 
cap. ly wi 
.06g2 syste 
per Hz 

1 of 4 






























































































N) 


Mtr, t.ilc 

Model 

U.llVIo 

Typo 

■ 

and 

Damping 

linearity 

Tliortiial 

Sen,*,, 

Sluit 

Thermal 

Dias 

ShlTL 

'/.oro 

Offset 

Trons- 
vorsc 
Sens' V 

Noise 

QLy 

Price 

Si VC or 
V.'eiglit 

Sundstrand 

■ 


Non- 

pendulous 

Servo 

1 v/g 

650 Hz 
±5% to 
390 Hz 

±0.05% FS 

1% 

per 

100°F 

0.05 g 

per 

100«F 

±0.1 g 

None to 
actual 
sen. 
axis 

5 mV 
RMS 

$595 

3.3 oz 

Ferranti 

FA. 21 

±20 g 

Pendulous 

Servo 

current 

S ma/g 
± 0.01% 

N.S. 

electronic 

or 

viscous 

±0.05% 

accuracy 

0.56% 

per 

100*F 


3xl0~'^g 




80 gm 

Sanders 
> 1000 hxs 

30-5 

±5 g 

Pendulous 
differen- 
tial xfmr. 

AC 

1.12Vrms 
per g 

33 Hz 
0.5 
fluid 

±0.2% FS 
to ±2.5 g 
±3% FS 
to ±5 g 

4% per 
100*F 


0.05% 

FS 




4.5 oz 

Setra 

113 

triaxia! 

±25 g 

Seismic 

disc 

cap. p.o. 

60 mV/g 

800 Hz 
0.7 gas 
squeeze 
film 

+1% FS 

±2% FS 
per 
lOO^F 

±2% FS 
per 
100“F 

±50 mV 
at 77“F 

O.Olg/g 

<2 mV 
rms 

1 - 4 
$142 
per axis 
C$425) 

4 oz 

Larson 

Cal 

±5 g 

Pendulous 

Servo 

1 V/g 

120 Hz 
0l4 to 1 
Electro 
magnetic 

±1% FS 
Lin, HYST, 
Repeat- 
ability 

±3% 

per 

100“F 



±1% FS 

0.005g/f 


Small 

$129 

1 oz 

Bell and 
Howel 1 

4-205- 

0001 

±1.0 g, 

±2 5 g, 
or ±5 g 

Unbonded 

strain 

gage 

40 mV. 
full 
range 
5 V ex. 

"High” 
0.7 ± 0.1 
@ 77“F 

±0.75% of 
full rang* 
lin. + 
hyst . 

temp. ,co 
-65®F to 

npensated 

+250®F 


"Low" 


Small 

$415 

4.5 oz 

West Coast 
Research 

703 

0.5, 1, 
2, 5, or 

10 g 

Bonded 

strain 

gage 

ImV/V TS, 
lOinV/V 
FS, IV 
FS, or 
5 V FS 

20 - 500 
Hz 
0.4 

% FS 
0.5 lln 
0.5 hyst 
0.25 re- 
peat. 

1% 

per 

lOO^F 

1% 
per 
100 “F 


1% 


1 

5295 to 
5340 
1000 
5177 to 
5204 

3 to 7 
oz 
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> 

I 

w 


MTr, Life Model Range Type 


•?i°xL%fr \»v 

axis units ^ 
avail. 


Bourns - 
10^ cycles 




Damping 


dry gas 
damping 


1 V/g 30 Hz 
10 V ex. 0.7+0. 3 


Gould 

Statham 


Schaevitz 


Donner 


^5 g 

Unbonded 

4 mV/V 

75 Hz 

±4 g 

strain 

Full 

0.7±0.1 


gage 

scale 

300 Hz 

+5 g 

Porce 

±5V±0.2% 

125 Hz 

Rebalance 

Full 

0.55 to 


Servo 

range 

0.75 

±20 g 

Pendulous 

±4 V 

ISO Hz 

to 

Servo 

Full 

1.2±0.2 

±40 g 


range 

@ 20®C 


Inertial 

grade 


Jnder 

develop- 

ment 


I Solid 
I State 


±10 V 

Full 

range 


700 Hz 
0.03 
squeeze 
film 


LINEAR ACCELEROMETERS 


Non^- , Tliermal Thermal Zero Trans— 
linearity Sens, Bias Offset verse 
Shift Shift Sers'v 


Size or 
Weight 
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EalLHIRS 


B 


FROM NATIONAL TECHNICAL INFORMATION SERVICE 



Development of Pre-Mining and Reclamation Plan Rationale 
for Surface Coal Mines 

PB-258 041/SET/PAT 590 p PC$22 50/MF$7.00 

Manual of Respiratory Protection Against Airborne ^ 
Radioactive Materials 

PB-258 0S2/PAT 147 p PC$6 00/MF$3 00 

Design and Construction of a Residential Solar Heating and 
Cooling System 

PB-237 042/ PAT 233 p PC$8.00/MF$3.00 
1976 Energy Fact Book 

ADA-029 331 / PAT 195p PC$1250/MF$1250 

Impacts of Construction Activities in Wetlands of the United 
Stdtos 

PB-256 674/ PAT 426 p PC$11 75/MFS3 00 

Standardized Development of Computer Software. Part I: 
Methods 

N76-30-849/ PAT 389 p PC$10 75/MF$3 00 

A Methodology for Producing Reliable Software, Volume I 
N76-29-945/ PAT 228 p PC$8 00/MF$3.C0 


Flovr and Gas Sampling Manual 

PB-258 080/ PAT 102 p PC$5 50/MF$3 00 

Solar Heating and Cooling in Buildings: Methods of Economic 
Evaluation 

COM-75-1 1070/ PAT 48 p PC$4 00/MF$3 00 

Data Base Directions The Next Steps 
PB-258 103/PAT 177 p PC$7 50/MF$3 00 

Comparative Study of Various Text Editors and Formatting 
Systems 

ADA-029 050/ PAT 93 p PC Not Available/MF$3 00 

Explaining Energy: A Manual of Non-Siyie for the Energy 

Outsider Who Wants in 

LBL-4458/PAT 78 p PC$4.50/MF$3 00 

A Survey of State Legislation Relating to Solar Energy 
PB-258 235/ PAT 166 p PC$6 75/MF$3 00 

Cost Estimating Handbook for Transfer, Shredding and 

Sanitary Landfilling of Solid Waste 

PB-256 444/ PAT 85 p PC$5 00/MF$3 00 

Coat Liquefaction Design Practices Manual 
PB-257 541 / PAT 372 p PC$1 0.50/M F$3 00 


HOW TO ORDER 


When you indicate the method of pay- 
ment, please note if a purchase order is not 
accompanied by payment, you will be billed 
an additional $5 00 ship and bill charge And 
please include the card expiration date when 
using American Express 

Norma! delivery time takes three to five 
weeks It is vital that you order by number 


or your order will be manually filled, insur- 
ing a delay You can opt for airmail delivery 
for S2 00 North American continent, S3 00 
outside North American continent charge per 
Item Just check the Airmail Service box If 
you’re really pressed for time, call the NTIS 
Rush Handling Service (703)557-4700 For a 
$10 00 charge per item, your order will be 
airmailed within 48 hours Or, you can pick 
up your order m the Washington Informa- 
tion Center & Bookstore or at our Springfield 
Operations Center within 24 hours for a 
$6 00 per item charge 


You may also place your order by tele- 
phone or if you have an NTIS Deposit Ac- 
count or an American Express card order 
through TELEX The order desk number is 
(703) 557-4650 and the TELEX number is 
89-9405. 

Thank you for your interest in NTIS. We 
appreciate your order 


METHOD OF PAYMENT 


Q Charge my NTIS deposit account no. 

n Purchase order no 

r~l Check enclosed for $ 

□ Bill me. Add $5.00 per order and sign below, (Not avail- 
able outside North American continent.) 

Q Charge to my American Express Card account number 


NAME 

ADDRESS 

CITY STATE ZiP. 


Card expiration dale 

Signature 

□ Airmail Services requested 

Clip and mad to 


National Technical Information Service 
U.S DEPARTMENT OF COMMERCE 
Springfield. Va 22161 
(703) SS7-4650 TELEX 89-940S 


Item Number 

Ouaniiiy 

Unit Price* 

Total Price’’ 

Paper Copy 
(PC) 

Microfiche 

(MF) 





















All prices subject to change The prices hub Total 

above are accurate as of 4/77 additional Charge 

Enter Grand Total 





Foreign Prices on Request 




A-4 


Mfr, & Life 


Range 

Type 


Damping 

Non- 

linearity 

Thermal 

Sens. 

Shift 

Thermal 

Bias 

Shift 

Zero 

Offset 

Trans- 
verse 
Sens ' V 

Noise 


Size or 
Weight 

Conrac 

24185 

i 

±1 g 

to 

±20 g ■ 

Spring 
restrained 
seismic 
mass, pot 

i 

±5 V 
Full 
range 
10 Vex. 

8 to 41 Hz 
0.3 min. 

+1% FS 
with 

0,2g 60Hz 
vibration 

±l^addi 
error t 
-65®F 1 

tional 

ver 

0 200*F 

1 to 2% 
@ 

25®C 

O.OOSg/fe 

1 

i 

1 

1 

i 


6 oz 
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Manufact'r 
and Life 

Model 

Range 

Type 


f and 
m 

Damping 


Thermal 
Sens . 
Shift 

Linear 

Acceler- 

ation 

Sens. 

Thres- 

hold 

X Axis 
Sens. 

Noise 

Qty, 

Price 

Size or j 
Weight 1 

Donner 

4590 

■ 

Liquid 

rotor 

sensor 

±5 V 
full 
range 

30 Hz @ 
1 r/s2 

0.1% of 

full 

range 

3% per 
100°F 

0.02% 

per 

g 

0.001% 

full 

range 

O.Olr/s^ 

per 

r/s2 

0.07% 

full 

range 

1-10 

$3000 

100-500 

$2500 

3 lbs 

Donner 

4591 

±10-200 

rad/sec2 

Liquid 

rotor 

sensor 

±2.5 V 

full 

range 

To 150 Hz 
0.6+0. 2 
@ 72 »F 

0.1% of 

full 

range 

3% per 
100®F 

H 

0.001% 

full 

range 

5 mV per 
rad/sec2 

lOmV 

rms 

1-10 

$3000 

100-500 

$2500 

8'. 5 oz 

Schaevitz 

ASBP-’S 

-1C 

-15 

5 r/s^ 
10 r/s2 
15 r/s2 

Torque 

rebalance 

Servo 

±5V±2% 

full 

range 

12 to 24 
Hz 

0.55 to 
0.75 

±0.1% 

full 

scale 

1% per 
100“F 

O.l-r/s^ 

per 

g 

< 0.1% 

zero 

offset 


5 mV 
rms 

1-2 

$664 

100-199 

$504 

3 oz 

Donner 

4577 

LOOr/s^ 

torque re- 
balance 
servo 

0.025Vdc 
per xlsr 






0.004 V/ 
r/s^ 

lOmV 

500 

$400 

5 oz 1 

\ 

X 
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r 

Manufact’r 
and Life 

Model 

Range 

Type 

dc Sens- 
itivity 

fp and 
Damping 

■ 

Thermal 

Sens. 

Shift 

Linear 

Acceler- 

ation 

Sens. 

Thres- 

hold 

X Axis 
Sens. 

Then. 
Bias 
Shift 
5r Nol 

Qty. 

Price 

le 

Size or 
Weight 

Donner 

MTBF 

> 60,000 hrs 

8160 

jKs9|H 

■■■■AH 

Liquid 
rotor 
; with 
integrator 

±5 Vdc 

full 

scale 

^nl 

0.007 Hz 
fn2 
30 Hz 
0.6+0. 2 

< 1% 
full 
scale 

< 5% 

per 

100“F 

* 

Note +2 
slope 
below 
0.007Hz 

<o:ooi% 

full 

scale 

0.005 

°/sec 

per 

®/sec 


1-20 

$1455 

100 

$1185 

500 

$995 

< 10 oz 

1 

Humphrey 
3 units cunt, 
op for 6 yrs. 

RT03- 

0108-1 


Electro 
fluidic 
1 axis 


50 Hz 

±2“/e6C 

5% 

-40®F to 
+170“F 




10% 

per 

lOO'^F 

1-5 

$935 

2.6" D 
0.4" L 

Humphrey 
3 units cont. 
op for 6 yrs . 

ATIO- 

0102-1 


Electro 
fluidic 
2 axis 




S% 

-40°F to 
+170‘’F 


2% 

on S] 

allable 
ecial I 

10% 

per 

1007 

1-5 

$1675 

2.6" D 
4" L 

Humphrey 
3 units cont. 
op for 6 yrs. 

RT02- 

0201-1 

Pitch 

teovs 

Roll 

tSdO^/s 

iTaw 

t60“/s 

Electro 
fluidic 
3 axis 

0 to 
+ 5 Vdc 

30 Hz to 
80 Hz 

. +1% 

5% 

-407 to 
+1707 


orde: 
$50 1 

add / 
o 100 \ 

10% 

.per 

1007 

1-5 

$3500 

3" X 

4.6" X 1 
3.7" 1 

BAC 9 XiCd • 

RG408 

Dart 

t20°/sec 

to 

tSOO^/se 

2 axis 









1-5 

$3500 

0.7" D 
2.2" L 

Honeywell 

GG2500 

LC02 

t4807sec 

Magneto- 
hydro- 
dynamic 
with spin 
motor 
2 axis 

ISmVrms 
per Vsec 
±2% 

100 Hz 

i 

0.1% FS 


0.05 

deg/sec/ 

g 

0.01 

deg/sec 

0.5% 

FS 

150 
mVrms 
® IkHz 
band 
width 

1 

$1000 

w/o 

elec- 

tronics 

0.7" D 
1.8" L 
70 gm 









' 
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r 

Manufact’r 
and Life 

Model 

Range 

Type 

dc Sens- 
itivity' 

f and 
m 

Damping 

Non- 

linearity 

Thermal 

Sens. 

Shift 

Linear 
Acceler- 
ation 
Sens . 

Thres- 

hold 

X Axis 
Sens . 

Noise 

Qty. 

Price 

Size or 
Weight 

Honeywell 
• cont . 

with 

electronic 
package 
for above 

II 

plus 

EG1030 

AB03 

±100'/se 

It 

87mVdc 
per */se( 
±2% 

-3 db 
0 70 Hz 

j 

0.1% FS 


0.05 

deg/sec/i 

0.01 

deg/se< 

0.5% 

FS 

ZmVrms 
5 null 
145 
nVrms 
@ FS 

. 

1 

$200 

* 

l.OS" X 
2.27" X 
1.9" 

1 

1 

1 


• ^ i 

# I 

ANGULAR. RATE SENSORS 6 lQOt 7 onloB only pago 2 oi 2 












/ 

Manufact’r 

Model 

if Axis 

Type 

Range 

Bj 

Scale 

Factor 

Accuracy 

Bias 

Accuracy 

Noise 

Temp. 

Stabil- 

ity 

(dc) 

Sensi- 

tivity 

Qty. 

Price 

Size or 
Weight 

Linearity 

" Develco 

9200C 

3 

Fluxgate 

±600 mg 

±1® 

±1% 

±1% FS 

1 Gamma 

p-p 

@ IHz 
b.w." 

0.4 mg 
/“C 

■ 

1-3 

$600 

to 

$700 

1000 

$360 

6 oz 

±0. 5,% 
FS 

Superconduct- 
ing Technolog) 

F203 

3 

■ 

±600 mg 

±1“ 

±1% 

±1% FS 

NS 

■ 

2.5F per 
600 mg 

II 

3.5 oz 

±0.5% 

FS 

Schonstedt 

SAM - 
73C 

3 

Fluxgate 

±600 rag 

±1“ 

1 

±1% 

±0.4% 

FS 

1 Gamma 

p-p 

@ IHz 
b.w. 

0.14 

mg/®C 

2.5 V 
per 600 
mg 

1-10 

$1650 

1000 

$600 

5 oz 

±0.4% 

FS 

Infinetics 

Friskem 

Mark-26 

1 

Fluxgate 

±600 mg 

!»s 

NS 

NS 

NS 

NS 

0.192 V 
per 600 
mg 

1-3 

$44 

1000- 

2000 

$34 

1" D 
2.5'* L 

Reason- j 
able ■ 

Spartan 

Electronics 

RCI 

1 

Fluxgate 

±1500 rag 

NS 

i 

±2% FS 

NS 

0.4 nT 

p-p 

T = 60sec 

NS 

1.2 V 
per 600 
rag 

1-6 

$1000 

1000 

$200 

1.75" X 
1.5" X 
0.75" 

1% FS 

Schonstedt 

RM62-? 

2 

* Sii 

ilar to 
Mr. Uptoi 

SAM 73 C 
^117/76 

per Telic 

>n 




100 

$550 

1 lb 
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iiilPRODUCIBILITY OF THB 
ORIGINAL PAGE.B POOR 


Mfr. & Life 

Model 

Range 

Type 


*Static 

Error 

Thermal 

Bias 

Shift 

Thermal 

Sens, 

Shift 

Error @ 
Sea Level 
±50“F^ 

Time 

Cinstanl 

Initial 

Zero 

Balance 

Long 

Term 

Stab. 

Qty 

Price 

Size or 
Weight 

Setra 

B 

0 to 20 
PSIA 

Quartz 
capsule 
cap. p.o. 

5 V FS 

±0.06% 

0.1% per 
100®F 

0.1% per 
100 “F 

±35 
feet , 

10 ms 
max 



1 

$975 

9 oz 

Hamiliton 

Standard 


0 to 20 
PSIA 

Vibrating 

cylinder 

Freq. 
prop, to 
pressure 

+0.016% 

0 . 38% pel 
compensat 


±13 

feet 

Frequ 

linea: 

;ncy is n 
■ with pr 

3n- 

sssure 

1 

$1380 

6.4 02 

Rosemount 

1241A 

^9 

metal 
capsule 
cap. p.o. 

BBS 

+0.25% 

reading 

+20' 

0.75% of 
+40 feet 
static ei 

reading 

inc. 

ror 

±47 

feet 

Out] 

witl 

lut is li 
1 altitud 

lear 

1 

$1155 

24 oz 

Senso-Metrics 

3 X 10^ 
cycles 

SP 

65 E 

3 to 15 
PSIA 

Bonded 

strain 

gage 

200 mV 
FS 
±2% 

+0.05% 

0.5% per 
100*F 

0.5% per 
100“F 

±100 

feet 


2% 

FS 

max 


1 

$900 

3 oz 

Gould 

Statham 

PA 824- 
15 

0 to 15 
PSIA 

Beam 

diaphragm 

deposited 

strain 

gage 

30 mV 
FS 

+0.15% 

0.5% per 
100“F 

0.5% per 
100°F 

+109 

feet 

0.05 

ms 

< 2% 
FS 

0,05% 

FS 

1 

$1265 

14 oz 

Bell and 
Howell 

CEC 

1000 

0 to 15 
PSIA 

Diaphragm 

deposited 

strain 

gain 

30 mV 
FS 

+0.25% 

0.5% per 
100“F 

0.5% per 
100®F 

±123 

feet 


±2% 

FS 


1 

$450 

5 oz 

ELH 

Electronics 

DV 

0 to 20 
PSIA 


5 V FS 

+0.3% 

1% per 
100°F 

0.5% per 
100°F 

±178 

feet 

0.05 ms 

Adjust 


1 

$420 

17 oz 

Bourns 

200 - 
438 - 
1002 

Zero 

to 

35,000' 

Diaphragm 

with 

potentio- 

meter 

10 V FS 

+1.5% FS 

voltage 

ratio 

2.5% errc 
per 1£ 

r 

0®F 

±825 

feet 

Outpi 

wii 

t is lin 
h altitw 

far 

le 

1-10 

$375 

13 oz 
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L-10 


t/fi . & Life 

i\lociel 

Range 

Type 


*Static 

Error 

Tlicriiial 

Bias 

Shift 

Thermal 

Sens, 

Shift 

Error 
Sea Leve' 
±50'’FZff 


Initial 

Zero 

Balance 

Long 

Term 

Stab. 

Qty 

Price 

Size or 
Weight 

Schaevitz 

PTD - 
310A - 
300W 

0 to 18 
PSIA 

Diaphragm 
capsule 
with LVDT 

2.5V FS 

+1.1% 

Operatinj 
range -6J 
200®F 

temp. 
°F to 

±316 

feet 

1 ms 

Adjust 


1-2 

$542 

18 oz 

Tyco 

AB-15- 

ADE 

0 to IS 
PSIA 

Diaphragm 
with semi- 
C. strain 
gage 

100 mV 
FS 
±1% 

+0*5% 

1% per 
100“F 

1% per 
100®F 

±245 

feet 

, 




1 - 9 
$285 
100 
$142 

2 OZ 

Stmdstrand 

314 A 

0 to 15 
PSIA 

Servo 

0 to lOV 
FS 

.±0,057% 

0.4% per 
100*F 

0.5% per 
lOO^F 

±92 

feet 

10 ms 



1 

$810 

1.6" X 
1.6" X 
1.9" 

Viatran 

304 


Diaphragm 

with 

strain 

gage 

High 

level 

available 

±0.15% 





Adjust 




Consolidated 

Control^ 

415G40 


Diaphragm 

bonded 

strain 

gage 

30 mV 
FS 

±0.25% 

TherraaJ 
< 0.5% 
per 10( 

sens, 

FS 

op 






3.5 oz 

Gultoa 
"Long " 

3255 

0 to 15 
PSIA 

Aneroid 

capsule 

with 

LVDT 

0 to 5 V 
FS 

±0.5% 

Total eri 
ing temp, 
-65'F to 
±0.75% 

or includ- 
from 
2S0°F 

+ 175 
“ft. 



out- 

stand- 

ing 


5 oa 

Gulton 

3261 

0 to IS 
PSIA 

Aneroid 

capsule 

with 

pot 

0 to lOV 
FS 

±1% FS 

Thermal 
1% per 

sens. 

100°F 

±317 

feet 

0.3% 
repre: 
sea li 

pot reso 
ents 85' 
!vel 

lution 

at 


2.5 oz 

Daytronic 

502 

0 to 15 
PSIA 

Diaphragm 

with 

strain 

gage 

30 mV 
FS. 

±0.25% 

0.5% per 
100®F 

0.5% pel 
100“F 

±122 

feet 

0.3 ms 

±2% 


1 

$295 

2.5" D 
3.5" L 
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M.unU.icL 'r 
and Liio 

Mtniol 

ll.mtjc 

Type 


♦Static 

Error 

Theriiml 

Dins 

Shift 

Thermal 

Sens, 

Shift 

Error 0 
Sea Level 
iSO^FAT 

Time 

Constant 

Initial 

Zero 

Balance 


Qty 
; Price 

i 

Size or 
Weight 

Honeywell 


AIR 

mputer 

Silicon 

pressure 

capsule 

Gen 

CPU 

sral purpo! 
employed 

;e 


±14 

feet 

■ 

>r not 
>arately 

Alti- 
tude + 
air- 
speed 
+ roach 

1 

$15 k 

19 lb 

Conrat 
2 X 10^ FS 
cycles min. 

4715LC 


Silicon 
capsule 
and strain 
gage 

50 mV 
+2% FS 

B 

0.5% per 
100®F 


±131 
feet j 


±1% 

FS 

1 

1 

1.12" D 
1.62" L 

National 


0 to 15 
PSIA 

Solid 

state 

silicon 

+2.5 V 
to 12.5V 
FS 

+2.1% 

2.7% per 
100®F 


±800 

feet 


Inc. in 

static 

error 

■ 

1 

1 

1 

$80 

5 gm 

Cognition 

CG4208 
1 

0 to 30 
inches 
Hg 

Solid 

state 

silicon 

+1.5 V 
to +4.5V 
FS 

±0.2% 

1 


effect 

100"F 

±140 

feet 


0.2% 

0.3% 

g 

2.5" X 
1.5" X 
0.75" 

Kavlico 

■ 

0 to 15 
PSIA 

26 Vac 
LVDT 
400 Hz 

2 Vac 
+3% FS 

' 

1% per 
lOO^F 

1% per 
lOO^F 





■ 

2.25 L 

1.25 D 

Transducer 
Systems Inc. 

ac/ac 

or 

dc/dc 

0 to 15 
PSIA 

LVDT 1 

j 

100 mV 
rras or 
2 Vdc 
FS 

+2% 

1% per 
100“F 


±1400 

feet 


j 


1 

$150 ac 
$250 dc 



Hr; 

D 

IV 

2" 

1 


3V 

Robinson 

Halpem 

155 

0 to IS 
PSIA 

Capsule + : 
LVDT + i 
electronic 
circuit 

0.5 to 
5.5 Vdc 
FS 

+0.5% 

FS 

Compensal 
±2% FS oi 
cal 0" t( 

ed to 
[ 75'F 
1 180°F 

+316 

feet 


Adj, on 
zero and 
span 


1 

$175 

8 oz 

Rosemount 

1332A 

0 to 15 
PSIA 

Capsule 
cap. p.o. 
el ectxonicj 

0 to 5V 
FS 

±0.11% 

Combined 
effect 
per 100°F 

thermal 

±0.5% 

±77 

feet 

5 ms 
max 

’ 1 
i 


1 

$445 

7 oz 
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Initial 

Zero 

Balance 


-45“F 

to 

+160®F 


+10mV 




QLy 

S 1 y.a or 

I’l Ice 

V/ci^ht 

1 to 5 
$2360 

2 lb 

12-24 


$2125 



1-9 

135. 
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A-13 


Viaimf act' 1 
and Life 

Model 

Kange 
PS ID 
(knots ) 

Type 


Static 
Error 
PS ID 

Xliermnl 

Bias 

Shift 

Thermal 

Sens, 

Shift 

I: 

Kf! j wksI 

Initial 

zero 

Balance 

Over 

Pres- 

sure 


Qty 

Price 

Siao or 
h'oi(;ht 

Sensotec 

TJE 

2 

(285) 

Bonded 

Strain 

Gage 

20 mV 
FS 

+0.1% 

0.25% 

per 

100®F 

0,25% 

per 

lOO^F 

3 

(342) 



■ 

Bj| 

2.5" L 
2.25" D 

Sensotec 

Z 


Bonded 

Strain 

Gage 

20 mV 
FS 

+0.3% 

0.5% 

per 

lOO'F 

0.5% 

per 

100°F 

3 

(342) 




1 

$505 

2.5" L 
2.25" D 

Sensotec 

AS 

2 

(285) 

Bonded 

Strain 

Gage 

20 mV 
FS. 

±0.5% 

1%" 

per 

100“F 

1% 

per 

lOO'F 

3 

(342) 




1 

$410 

2.5" L 
2.25" D 

Sensotec 

AlO 

2 

(285) 

Bonded 

Strain 

Gage 

30 mV 
FS 

+1% 

1.5% 

per 

100°F 

2% 

per 

lOO'F 

3 

(342) 




g| 

2.5" L 
2.25" D 

Celesco 

P70 


Diaphragm 

Variable 

Reluctance 

25 mV/V 
at 

3000 Hz 

±0.7% 

o£ 

^nge 

< 1.7% 
error 

■1C 

Thermal 

per 

0®F 



200% 

Of 

range 



14 oz 

C.J. 

Enterprises 

CJVR 

±0.1 

to 

±S00 

Diaphragm 

Variable 

Reluctance 

so mV/V 
at 

5000 Hz 

+0.7% 

of 

Range 

1.8% 

error 

1C 

Thermal 

per 

0°F 



200% 

of 

range 



15 oz 

Gould 

Stratham 

PL 872 
-S 

91 

Diaphragm 

Deposited 

Strain 

lage 

30 mV 
FS 

+0.2% 

0.5% 
per 
100 ®F 

0,5% 

per 

100®F 


+2% 

FS 

200% 

of 

range 


ll^llllllll^ 

14oz 

Bourns 


(40) to 
(400) 

Diaphragm 

with 

Potentio- 

meter 

10 V 
FS 

+2% 

FS 

Voltage 

Ratio 

1.5% e 
per 1 

rror 
00 “F 

Outpi 

wit 

t is lin 
h alrspe 

ear 

ed 



13 oz 

Schaevitz 

PTD - 
310 D - 
50 W 

1.84 

(270) 

Diaphragm 

Dapsule 

with 

LVDT 

2.5 V 
FS . 

±1.1% 

Operat 

temp. 

-65“F 

ing 
range 
to 200+F 


Adjust 



1-2 

$445 

18 oz 
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> 

I 

!-> 


J.ljiuiract ’ z' 
and Life 



Diaphragm 

with 

Strain 

Gage 


Consolidated 41GB67 
Controls 



Inconel-X 

Diaphragm 

with 

LVDT 



0 to 5 
(435) 

Aneroid 

capsule 

with 

LVDT 

5 

(435) 

Int . ckt 
silicon 
diaphragm 
diffused 
strain gage 


5 

3S) 

Diaphragm 

with 

Strain 

Gage 
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M.umfact. 'x 
and Life 

Aiodcl 

Uangc 

PSID 

Otnotrf) 

Type 


* 

Static 

Error 

PSID 

Thermal 

Bias 

Shift 

Thermal 

Sens, 

Shift 

KinHSXI'li 

Initial 

zero 

Balance 

Over 

Pres- 

sure 

■ 

Qty 

Price 

Size or 
Weight 

, Conrac 

4715H- 

D 

10 

(1400) 

Silicon 
capsule 
and strain 
gage 

5 V 
FS 

±0.4% 

0.5% 

0.5% 


±1% 

FS 

2 X 

range 

1 

1 

$1500 

35 

$700 


Nation.al 

LX1501D 

-5 to 5 
(435) 

Solid 

state 

silicon 

+2.5 V 
to 

+12.5 V 
FS 

6% 

of 

5 PSID 

6% 

per 

100“F 

6% 

per 

100“F 


inc. in 

static 

error 

40 

PSID 

■ 

1 

$85 

5 gra 

Transducer 
Systems Inc. 

ac/ac 

dc/dc 


LVDT 

100 mV 
RMS or 
2 Vdc 
FS 

±2% 

1% 

per 

lOO^F 





-25®F 

to 

300®? 

op. 

1 

41 Cn ar 

|ac fdc 

nffiVMI 




Robinson- 

Halpem 

150 

3 

(342) 

Capsule + 
LVDT + 
electronic 
circuit 

0 to 1 V 

±1% 

Cal. 
and li 
in 

ca] 

at 30°F 
0"F with- 
% of 75°F 


Factory 
set 
±1% FS 

1.5 X 
range 

-40“F 

to 

220“F 

op. 

1 

$130 

8 oz 

Rosemount 

1221B3 

2.456 

(300) 

Capsule 
cap p.o. 
elec- 
tronics 

6.14 V 
FS 

±0.1% 
inc , in 
op. • accy. 

operati 
accurac 
-5S-C t 

ng 

y ±0.5% 
0 +71*0 


Inc. 
in op. 
accuracy 



1 

$785 

13 oz 

Rosemount 

S42K2 

High 

level 

dc 

outputs 

lAltitude: 

(airspeed; 

[dynamic 

h 

IAS 

pressure: 

-1000‘ 

75 

4C 

’To +40, C 
to 340 1< 
0,27 to 

DO* 

nots 

6.29" Hg. 

±(40* + .7 
±(2.5K+ . 
±(.S% FS- 

X readin] 
4% readii 
+ 1% read: 

!) )-45“i 
ig) ? to 
Lnsy+IBO" 

.2mV/ 

ft 

P- 20mV 
K 

1V/"H 

1-5 
$2360 
> 12-24 
$2125 

y 

9 

2 lb 
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'sOOi II 30V5 

mi m xa^iTiaionaosda^: 


Maiiulact' 1 
and Life 

Model 

1 Range 

Type 


Non- 

linearity 

Thermal 

Bias 

Shift 

Thermal 

Sens, 

Shift 

■ 

* 

Static 

Error 

Band 



Qty" 

Price 

Size or 
Weight 

Celesco 

TCC - 
PT - 
101 

Up to 
500*' 

1 

Cable pot 
spring ' 

return 

10 V FS 

< 0.1% 








5V’ X 

2" X 
25s" 

Bouttjs , 
40 X 10° 
cycles 

194 

i 8" 

i 

1 

! 

Linear 

pot 

10 V 
FS 

±0.5% 1 

1 

Temp 

-65-F 

Range 
to ZSQop 

1 lb 
max 

1 

■ 



Hi 

0.75" D 
11" L 

Bourns - 
40 X 10° 
cycles 

i-163 


Linear 

pot 

10 V 
FS 

±0.5% 

Temp 

-es^F 

Range 
to 350*F 

1 lb 
max 

i 

i 


! 

Sealed 

1 

$255 

25 

$135 

0.375" D 
By 3,4,5 
or 6" 

BLH 

Electronics 
10^ cycles 

416347 

4i6348 

416349 

423977 

.75 and 
2" 

Resistive 

22.5 to 
45 mV 
FS 


2% 

per 

100®F 

0.5% 

per 

lOO'F 

25 gm 

0,25% 
CO. 5% 
for 2") 


+15®F 

to 

11S“F 

op. 

range 

1 

$145 
.1-.75" 
$225 2" 

8-12 oz 

West Coast 
Research 

557 

±0.01" 

to 

410" 

Beam 

deflection 
with strai^ 
gage 

3 raV/V 
to 5 V 
: FS 

0.05% 

Thei 

Sensit: 

0.2% 

mal 

vity 

per 106®F 






1" 

cube 

West Coast 
Research 

557 L 

j 

0 to 
0.5" or 
+0.25" 

1 

same 

i 

20 mV 
FS 

0.75% 

0.5% 

per 

lOO^F 



2% 

i 

“40 to 
350®F 
expo- 
sure 
range 


1" X 
1,5" X 
3.5" 
3/8" 
shaft 

All tech 

1 

1800 

Up to 
300" 

Cable pot 

spring 

return 






Up to 
0.05% 





Kavlico 
50,000 Hrs 

Many 

±0.05" 

to 

6" 

stroke 

LVDT 

400 Hz to 
3 kHz 
IV/in/V 
typ. 

0.4% 

typ. 

1 


0.1% 

per 

100®F 

4 oz 
typ. 


5“ phase 
30 mV 
null typ 

-420®F 

to 

M50“F 
typ. op 

♦ 

Many 
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Mjuitroct'i 
nnd Life 


Humphrey 
10^ cycles 


The Accro 
Co. 


Hewlett 

Packard 


Transducer 

Systems 

Inc. 



RP93- 

0101-1 


Series 

1850 



7DCDT 

and 

24DCDT 




llange 

Type 

0.25" 

to 

1" 

stroke 

Linear 

pot 

2" 

to 

180" 

cable pot 

spring 

return 

12" to 
60" . 

Ultrasonic 

Magnetic 

+0.25" 

to 

±1" 

LVDT 

±0.05" 

to 

±3" 

LVDT 

+0.05" 

to 

+3" 

LVDT 

+1" to 
+6" 

others 
down to 
±0.005" 

LVDT 



Non- Tliernial 

linearity 

Shift 


Size or 
Weight 




0.044 

to 

460mV/in 


10 V FS 
or 

digital 


AC 

1.6 to 

4.6 mV/V 
DC 0.4 

mV/V 


4.3 to 
100 V/in 


±0.4 V +0.25% 

or or 

±10 V FS +0.1% 


1 to 2% 
others 
0 . 1 % 


1.2 to 2" 
10 to 16" 


23 to 

220 gm 


0.75" D 
3 to 14" 
Length 


LIHMR DXSPIACEHBNT TRANSDUCERS 


nonlinearity, hysteresis, nonrepeatability 
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MiimiXa( t’l' 
and L\Cc 

Model 

llange 

Tyj )0 

■■■■■ 

Non- 

linearity 

Tliermnl 

Dias 

Shift 

Tlioniial 

Sens, 

Shift 

Actuation 

Force 

*Static 

Error 

Band 



Qiy 

Price 

Si /.a or 
Weight 

Schaevitz 

Many 

±0.05 *' 
to 
+10" 

LVDT 

ae 0.1 
to 6inv/\ 
per mil] 

0.05% 

to 

4% 



contact- 
ing free 
slug 



h~65°F 
to 
300 “F 
op. 

1 

$32 to 
430 
6-11 
$28 to 
387 

4 to 580 
gm 

1 to 30 
inches 

Robins'on 

Halpem 

220 

Series 

+0.005" 

to 

+18" 

LVDT 

mm 

+0.07% 

to 

+0.5% 

0.3% to 
10% per 
100“F 


'Ion 

contact- 
ing free 
slug 


60 Hz to 
20 kHz 
units 



3/S or 
7/8" D 
0.5 to 
18" L 

Aviation 

Electric 

Ltd, 

1303 

1343 

1353 

1315 

1362 

1.27imn 

2.5 

10 

25.4 

76.2 

LVDT 

+5 V 
+10' V 
±5 V 

II 

+3 V 

+0.15% 

to 

±1% 



Spring 
loaded or 
free slug 
or self- 
aligning 
hearing 


Linear p 
also ava 
12.5 to 
mra Strok 

)tS 

Liable 

>540 

> 
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APPENDIX B. 


ADDITIONAL ANGULAR ACCELERATION DISCUSSION (Ch. 4) 

The angular accelerometers available measure the torque to accelerate 
an inertia of rigid or liquid matter. In principle, linear accelerometers 
can be used differentially. If only two are used, they are sensitive to 
an angular velocity vector component in the plane of their separation and 
sensitive axes, A set of four at right angles can be used to cancel this 
effect, or three can be deployed in a plane with the input axes 120 degrees 
apart . 

The measurements for two units as shown in Fig. 13-1 are 



= {f . g] 

y'^ = {f-b+[(wXR+wx (WxR)] • 6} 

For the ideal case with a and b parallel and each perpendicular to 

R and K = K =1* 
a b 

Id 

y + y = Rq ’2 + ^“ 1^3 • 

Thus the signal that is the sum of the accelerometer outputs measures 
angular acceleration but has error terms proportional to the products of 

B-i 



the two angular velocities. If four are used as shown in Fig. B-2, there 
is no error. 





Fig. B-2 


Similarly, for three units in a plane as shown in Fig. B-3 



Figure B-3 


with 


B-2 




each measurement xs of the form 

= (f + + (A) X w X R^) • a . 

This confxguration rejects the effect of linear specific force as before. 

Thus, three, four or more accelerometers can be located in a plane 
in a configuration that is insensitive, in principle, to angular velocity, 
but two accelerometers are not adequate due to their intrinsic sensitivity 
to angular velocity. 



F^r 10 millirad angle information accuracy, the angular accelerometer 

performance requirement depends upon the amount of time that errors are 

permitted to grow. Rapid' roll angular determination with a rolling mode 

time constant of a few lOths of a second would correspond to perhaps 1 sec 

smoothing time, whereas spiral and phugoid modes may require 100 times as 

_2 

long, thus the minimum error or undetected rate would be from 2 X 10 

- 6,2 , 

to 2 X 10 rad/sec . Assuming a maximum rolling rate of about 180 deg/sec 

and a rolling mode time constant of 0.1 sec, the maximum angular accel- 

' , 2 

eration would be approximately 30 rad/sec . We can calculate a dynamic 

,2 ■ 

range for the instrument as 0.01 to 10 rad/sec for an absolute minimum, 

. -3,2 -6 

0.4 X 10 to 30 rad /sec would be more desirable, and 10 to 100 

would give more information than is needed. Assuming a 10 cm separation 

between accelerometers, the equivalent linear acceleration uncertainty 

-4 

for even the minimum range is 10 g and this may be below the level of 
some accelerometers that would be acceptable for measuring the linear 


B-3 



In addition to basic sensitivity, the alignment, scale factor match- 

“6 

ing, and linearity must be in the 10 range to reject the 1 g specific 
force. 

Angular accelerometers, therefore, should be special purpose devices 
rather than a set of linear accelerometers for the aircraft autopilot 
measurement application. 
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